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SUMMARY 

A  method  for  predicting  the  propulsion-induced  aerodynamics  of  a  VSTOL 
aircraft  in  the  transition  flight  regime  has  been  developed.  This  method 
represents  the  formulation  of  wind  tunnel  test  data  consisting  of  pressure 
coefficients  at  numerous  locations  on  a  flat  plate  due  to  a  circular  jet 
exhausting  normally  into  crossflows  of  freestream  to  jet  velocity  ratios  from 
0.1  to  0.3.  As  a  result,  the  method  is  applicable  to  low-wing,  circular  jet 
subsonic  V/STOL  configurations  with  normally  exhausting  jets. 

This  report  documents  the  development  of  the  formulation  along  with 
validation  results  for  various  V/STOL  configurations  indicating  its  applica¬ 
bility  as  well  as  limitations.  The  computer  code  used  in  the  validation  of 
the  method  is  also  presented  along  with  the  required  input  and  configuration 
modelling  procedures. 
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INTRODUCTION 

The  V/STOL  transition  flight  regime  covers  the  velocity  range  from  low 
speed  hover  to  the  sustained  velocity  required  for  fully  wingborne  flight.  Jet 
exhaust  effects  have  a  major  influence  on  aerodynamic  characteristics  in  this 
regime  resulting  in  large  forces  and  moments  being  induced  on  the  aircraft. 

These  jet  effects  involve  the  interaction  of  the  jet  and  free  stream  flow  and 
are  primarily  the  result  of  four  viscous  flow  phenomena:  jet  blockage,  wake 
separation,  entrainment,  and  vortex  generation.  Due  to  complex  viscous  inter¬ 
action  of  these  four  phenomena,  the  individual  effect  of  each  component  has  not 
been  separated,  making  analytical/numerical  prediction  most  difficult. 

Current  prediction  techniques  applicable  to  the  transition  flight  regime 
consist  of  potential  flow  analysis  computer  programs  with  empirically  derived 
adjustments  or  corrections  to  account  for  the  viscous  interactions  which  char¬ 
acterize  transition  aerodynamics.  These  techniques  are  extremely  time  consuming 
with  respect  to  modelling  the  aircraft  configuration  with  panels  required  by 
the  potential  flow  analysis,  and  also  with  respect  to  the  computer  time  required 
to  run  the  program.  Unfortunately,  this  time  expended  does  not  necessarily 
result  in  an  associated  increase  in  accuracy.  As  a  result,  these  techniques  are 
primarlv  applicable  to  detailed  design  type  of  analysis,  being  too  expensive  for 
use  as  a  preliminary  design  tool. 

The  methodology  presented  here  represents  the  development  of  a  prediction 
technique  which  is  designed  for  use  as  a  preliminary  design  tool.  The  induced 
lift  and  pitching  moment  are  calculated  requiring  only  the  configuration  geom¬ 
etry  and  jet  location.  Since  most  of  the  available  experimental  data  applicable 
to  the  transition  flight  regime  resulted  from  specific  aircraft  configurations, 
it  was  decided  to  use  the  more  fundamental  data  obtained  from  tests  of  jets 
issuing  from  flat  plates.  This  data  formed  the  basis  of  the  prediction  tech¬ 
nique  which  could  then  be  applied  or  adapted  to  most  aircraft  configurations. 

Prior  to  the  formal  publication  of  this  report,  this  method  was  used  as 
the  basis  for  a  more  complete  prediction  technique  for  V/STOL  transition  aero¬ 
dynamics  with  no  restrictions  concerning  aircraft  configuration.  Coefficients 
are  applied  to  the  basic  lift  loss  value  calculated  by  the  method  contained 
herein  to  account  for  various  configuration  effects  including  wing  aspect  ratio, 
longitudinal  position  of  the  jet,  wing  height,  lateral  spacing  of  multiple  jets, 
nozzle  configuration,  and  jet  deflection  angle.  A  complete  discussion  of  the 
development  of  the  more  complete  technique  is  contained  in  reference  (6). 
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METHOD  DEVELOPMENT 

The  preliminary  design  stage  of  a  V/STOL  aircraft  involves  the  assessment 
of  various  conceptual  designs  to  determine  which  concept  most  satisfies  the 
design  requirements  to  justify  further  development.  As  part  of  this  assessment, 
the  propulsion  induced  aerodynamic  effects  must  be  predicted  to  ensure  the 
selection  of  that  design  which,  when  weighed  against  the  other  design  require¬ 
ments,  minimizes  the  negative  effects  and  maximizes  the  positive  effects. 

In  developing  a  method  to  predict  the  propulsion  induced  effects  in  such  a 
design  environment,  the  design  information  required  by  the  method  must  be 
limited  to  that  which  would  be  known  at  that  stage.  Accordingly,  the  present 
method  was  developed  under  the  guidelines  of  requiring  only  general  configura¬ 
tion  geometry,  jet  location,  and  free  stream-to-jet  velocity  ratio  for  the 
desired  flight  conditions.  The  approach  then  taken  to  develop  this  generalized, 
preliminary  design  stage  type  of  empirical  formulation  to  predict  the  transition 
aerodynamics  was  to: 

a.  Obtain  available  pressure  coefficient  data  measured  on  a  large 
flat  plate  due  to  a  jet  exhausting  perpendicular  to  the  free  stream. 

b.  Derive  an  expression  for  the  pressure  coefficient  data  as  a 
function  of  velocity  ratio,  V  ,  and  longitudinal  and  lateral  distances  from  the 
jet,  X/D  and  Y/D  respectively. 

c.  Use  the  finalized  expression  as  a  basis  to  integrate  over  the 
configuration  planform  area  of  interest  relative  to  nozzle  location  to  obtain 
the  jet  induced  aerodynamics  in  transition  flight. 

d.  Validate  the  resulting  method  initially  using  simplified  flat  plate 
or  low  wing  configurations  followed  by  more  realistic  complex  designs. 

e.  Modify  or  adjust  the  expression,  based  on  results  of  step  (d)  to 
account  for  configuration  variables  such  as  jet  nozzle  location  relative  to  wing 
and  fuselage,  planf orm-to-jet  area  ratio,  jet  deflection,  and  fuselage  contour. 

The  pressure  coefficient  data  used  to  develop  the  empirical  formulation  was 
obtained  from  flat  plate  data  generated  by  Fearn,  reference  (1).  This  data  was 
obtained  for  free  streara-to-jet  velocity  ratios  ranging  from  0.1  to  0.45  using  a 
four-inch- diameter  jet  exhausting  normally  from  a  flat  plate  24  jet  diameters 
wide  by  27  diameters  long.  The  pressure  coefficient  was  defined  as 

Cp  =  E_~  JjfL- 

<u 

where  p  is  the  difference  between  the  pressure  measured  at  the  test  condition  and 
the  pressure  measured  with  the  power  off,  and  p^  is  the  static  pressure  of  the 
free  stream  fluid.  Thus,  the  data  represents  the  induced  effects  produced  by  the 
jet.  As  plotted  in  figures  (1)  to  (4)  the  data  approximates  an  exponential 
trend  with  X/D  which  tends  to  flatten  out  as  Y/D  is  increased,  with  a  non¬ 
linear  peak  value  variation  with  Ve.  In  developing  an  expression  to  fit  this 
data,  each  of  these  variations  must  be  accounted  for  plus  the  variation  of  the 
X/D  location  of  the  peak  values  for  a  particular  Vg,  as  emphasized  in  Figure  (4). 
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The  equations  representing  the  induced  pressure  coefficient  data  were 
developed  in  essentially  two  parts.  The  first  part,  Cpu^,  equation  (1),  cal¬ 
culates  the  peak  pressure  coefficient  as  a  function  of  Y/T^  and  Ve.  The  second 
part,  equation  (2),  represents  the  remaining  normalized  data  points 

(normalized  by  the  peak  Cp  value,  Cp^^,  for  each  Ve).  The  two  equations 
are  then  simply  multiplied  according  to  equation  (3)  to  obtain  the  desired 
induced  pressure  coefficient. 


CpMAX 


_ -4.25 _ 

2 

(4V  -  1)  (3.25  V  +  1.4) 

e  (Y/T>e+  .5) 


(1) 


Cp 

NORM 


1 


[K,  (X/D-F)]  K2 
e  e 


(3.67  Y/De+  5)  Vg  (2) 

e  (X/De  +  4  Y/De+  2. 5)2 


F  =  (2.48  y/De-  1 . 6 )Ve  -  .1  Y/De-  .07 
and  for  (X/De-  F)  <_  0 

K1  =  _ 2 _ 

(-2.28  Ve  +  1.36)  Y/Dg 

K2  =  16  Ve  +  1.55  In  (Y/De)  -  1 

Y/De 

whereas  for  (X/Dfi-  F)  0 

K1  *  _ I _ 

1.1  Y/De 

K-,  =  -.13  (Y/De-  3 . 5  )2  +  1.8 

C _  =  CD  Cd  1,2  ^ 

p  PMAX  rN0RM 
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In  equation  (1),  the  constant  -4.25  establishes  the  maximum  value  which 
attains.  The  term  e(4Ve-l'1‘-  accounts  for  the  exponential  variation  of  Cpy,^ 
values  as  shown  in  Figure  (2),  with  the  term 

(Y/De+  .5)  (3-25  Ve  +  l-4  > 

providing  the  variation  in  values  with  Y/7^,  and  Vg  which  tends  to  flatten 

the  Cp  curves  as  Y/I^  increases,  as  seen  in  Figures  (1)  to  (4). 

The  first  part  of  the  Cp\j0R!g  equation 


_ 1 _ 

elKi  (X/Pe-  F)]  K2 

establishes  the  basic  exponential  shape  of  the  curve,  being  a  version  of  the 
general  equation  1  Replacing  the  x  term  by  CX/Be  -  F)  accounts  for  the 

e(ax)b 


X/I^  shift  in  the  CpJi^AX  values  with  varying  Vg  and  Y/D^  as  indicated  in  Figures 

(1)  to  (4).  The  term  K]_  represents  the  coefficient  of  x  required  to  curve  fit 
the  data  and  provides  a  corrective  effect  to  insure  the  exponent  remains  posi¬ 
tive.  The  overall  exponential  shape  of  the  curve  is  made  svmmetrical  about 
the  ordinate  axis  by  the  factor  Kt  which  again  varies  with  Vg  and  Y / Dg. .  The 
second  term  of  the  cp^gp^  equation 

(3.67  Y/I^+  5)  V4 

e(X/D  +  .4  Y/T)  +  2. 5) 2 
c  e  e 

modifies  the  original  shape  of  the  curve  at  high  negative  X/l^ values  to  account 
for  the  positive  pressure  coefficient  in  front  of  the  jet. 

With  the  equation  developed  in  this  form,  the  calculation  of  induced 
pressure  coefficients  is  relatively  simple.  However,  calculating  the  induced 
forces  for  an  entire  configuration  is  more  amenable  to  use  of  the  computer.  As 
a  result,  the  formulation  was  computerized  with  a  configuration  modelling  tech¬ 
nique  and  an  integration  procedure  included  as  part  of  the  code.  The  modelling 
technique  and  integration  procedure  are  discussed  in  Appendix  A.  A  listing  of 
the  computer  program,  a  discussion  of  the  required  input,  and  a  sample  output 
are  then  presented  in  Appendix  R. 
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METHOD  VALIDATION 

Since  the  data  used  for  formulation  represents  pressures  on  a  flat  plate 
due  to  normally  exhausting  jets,  the  method  mainly  applies  to  low  wing  config¬ 
urations  with  a  normally  exhausting  jet.  The  jet  should  be  centrally  located 
in  the  fuselage  or  wing  undersurface,  away  from  wing  or  fuselage  edges  which 
would  result  in  additional  circulation  effects  not  contained  in  the  data. 
Additionally,  nozzle  pressure  ratios  iuring  the  test  indicate  the  method  to  be 
applicable  to  subsonic  configurations.  Therefore,  the  method  was  validated 
against  configurations  commensurate  with  its  capability  plus  additional  con¬ 
figurations  to  indicate  possibilities  for  extending  its  capability. 

The  first  correlation  of  the  method  was  done  with  data  from  reference  (2) 
for  a  rectangular  wing.  Figure  (5)  shows  fairly  good  agreement  for  both  lift 
and  resulting  pitching  moment,  especially  since  the  data  contains  some  circu¬ 
lation  effects  causing  the  apparent  over-prediction  of  induced  lift  loss  while 
the  method  does  not  contain  circulation  effects. 

Figure  (6a)  shows  the  correlation  of  data,  also  from  reference  (2),  for  an 
elongated  body  with  a  jet  located  40%  from  the  leading  edge.  Excellent  results 
are  shown  for  jet  induced  lift  while  the  pitching  moment  exhibits  fair  correla¬ 
tion.  Additional  comparisons  with  the  same  elongated  body  but  with  the  jet 
located  60%  from  the  leading  edge  are  shown  in  Figure  (6b).  Very  good  agreement 
again  is  shown  between  the  test  data  and  predicted  results. 

For  a  delta  wing  configuration  with  a  small  jet,  also  from  reference  (2), 
excellent  agreement  with  test  data  is  shown  in  Figure  (7a)  for  a  velocity  ratio 
range  of  0.1  to  0.3.  This  velocity  range  is  representative  of  the  lower  and 
middle  transition  region  and  thus  a  good  correlation  in  this  range  is  considered 
to  be  of  primary  importance.  The  corresponding  pitching  moment  plot,  however, 
indicates  gross  disagreement.  This  is  due  to  questionable  data  since  the 
aberrant  points  on  the  plot  represent  a  large  hump  in  the  data  which  substan¬ 
tially  deviates  from  the  trend. 

Results  of  replacing  the  small  jet  with  a  larger  jet  in  the  same  delta 
wing  configuration  are  shown  in  Figure  (7b).  For  this  configuration,  the  pre¬ 
dicted  results  for  jet  induced  lift  are  slightly  greater  than  actual  test  data, 
but  again  this  difference  is  small.  Excellent  agreement  for  pitching  moments 
is  shown  between  the  predicted  results  and  test  data  for  the  velocity  range 
0.1  to  0.3. 

Further  validation  was  conducted  with  a  delta  wing  design  from  reference 
(3).  Figure  (8a)  shows  the  comparison  to  be  fair  for  the  configuration  with  a 
jet  to  planform  area  ratio  of  0.006,  whereas  increasing  the  ratio  to  0.024  pro¬ 
duced  quite  good  results  as  indicated  in  Figure  (8b). 

To  determine  the  possibilities  of  extending  the  method  application  to 
multiple  jet  conf igurations,  data  from  reference  (4)  for  a  four  jet  arrangement 
of  the  same  configuration  as  above  was  used  for  comparison.  The  method  was 
applied  in  the  same  manner  as  before  for  each  jet,  however  a  thrust  weighted 
summation  was  used  to  obtain  the  total  induced  force  as  given  by 

AL  =  / AL\  Ti  + 

T  W*  Ttotal 


’  AL 


4+1 


‘  ^  "'•'total 


(*k)  Tn 
/n  Ttotal 
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where  (AL/T)^  and  are  the  induced  lift  and  thrust  associated  with  an  individ¬ 
ual  jet,  and  TCota^  is  the  total  thrust  of  the  configuration.  Figure  (8c)  con¬ 
tains  a  planform  view  of  the  configuration  along  with  results  of  the  comparison, 
indicating  excellent  agreement  up  to  velocity  ratio  of  0.15,  where  the  predic¬ 
tion  begins  to  diverge  into  only  fair  agreement  at  Ve  =  0.25. 

A  final  validation  of  the  method  was  done  using  data  from  reference  (5) 
for  a  high  wing  configuration.  This  represents  an  additional  extension  of  the 
method's  applicability  beyond  the  low  wing  type  of  configuration.  As  shown  in 
Figure  (9),  the  comparison  is  quite  good  for  the  single  jet  configuration 
indicating  wing  height  to  be  a  secondary  effect  in  this  design.  However,  Figures 
(10a)  and  (10b)  indicate  an  important  limitation  of  this  method  being  the 
inability  to  account  for  additional  circulation  created  by  a  jet  located  near 
the  trailing  edge  of  a  wing.  The  additional  lift  and  subsequent  reduction  in 
induced  lift  loss  associated  with  such  a  jet/wing  arrangement  is  clearly  illus¬ 
trated  in  the  figure. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  method  developed  herein  to  predict  the  propulsion  induced  lift  of  a 
V/STOL  aircraft  in  the  transition  flight  regime  has  been  shown  to  be  quite 
effective.  When  applied  to  those  configurations  which  do  not  exceed  the  limi¬ 
tations  imposed  by  its  development,  the  method's  predictions  are  well  within 
the  accuracy  required  of  a  preliminary  design  phase  analysis  technique. 

However,  these  limitations  require  that  further  development  of  the  method 
be  directed  towards  extending  its  capability  to  include  other  configuration 
variables.  These  variables  include  nozzle  location  relative  to  the  wing,  jet 
deflection  angle,  and  sideslip  conditions.  Additionally,  nozzle  pressure  ratio 
effects  must  also  be  included  to  allow  prediction  of  supersonic  V/STOL  config¬ 
urations  which  are  typically  characterized  by  one  or  more  high  nozzle  pressure 
ratio  exhaust  jets. 
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Pressure  Coefficient  Distribution 


FIGURE  3.  Pressure  Coefficient  Distribution 
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FIGURE  5.  Comparison  of  Predicted  Results  with 
Test  Data  of  a  Rectangular  Wing 


FIGURE  6.  Continued 
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FIGURE  7.  Continued 
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FIGURE  10.  Continued 
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APPENDIX  A 


MODELLING  AND  CALCULATIONS  PROCEDURES 
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Modelling  the  Configuration 

A  simple  panelling  method  using  rectangular  segments  to  approximate  the 
geometry  of  the  planform  is  used  to  input  the  desired  configuration.  As  shown 
in  Figure  (A-l),  the  length  of  these  rectangular  segments  is  restricted  only  by 
the  size  of  the  configuration,  with  the  width  being  dictated  by  the  accuracy 
with  which  the  curved  and  angled  aircraft  components  are  to  be  modelled.  To 
further  refine  the  calculation  of  the  induced  pressures  within  each  rectangular 
segment,  an  integration  interval  is  used  and  required  as  input  by  the  computer 
program  which  defines  the  longitudinal  and  lateral  grid  of  rectangular  elements. 
The  mid-point  of  these  elements  defines  the  point  at  which  actual  pressure 
coefficient  calculations  are  made. 

The  only  caution  to  be  heeded  when  panelling  the  configuration  is  that  suf¬ 
ficiently  small  elements  are  defined  in  the  area  around  the  jet  to  enable  accurate 
calculation  of  the  rapidly  changing  induced  pressures  associated  with  this  area. 
Conversely,  modelling  of  the  planform  at  distances  beyond  five  jet  diameters 
from  the  jet  is  not  critical  since  induced  pressures  are  negligible  beyond 
this  distance. 

Once  the  configuration  is  panelled,  the  computer  program  is  structured  to 
calculate  the  induced  pressures  for  any  number  of  different  velocity  ratios. 

This  can  be  done  by  simply  inputting  the  number  of  cases  to  be  run  and  the 
associated  values  of  velocity  ratio;  a  more  detailed  discussion  of  which  is 
given  in  the  computer  program  section. 

Integration  Procedure 

The  integration  procedure  used  to  exercise  and  validate  the  developed 
method  consists  of  the  incremental  form  of  the  equation 

AL  =  2  £2 2  /  /  Cp  d(X/D)d(Y/D) 

T  if  qj  J  J 

The  procedure  used  involves  the  following  steps: 

(1)  Pressure  coefficients  calculated  at  various  longitudinal  stations  for 
a  constant  span  station  are  plotted,  with  the  area  under  the  curve  calculated 
assuming  the  pressure  coefficient  to  be  constant  across  the  various  longitu¬ 
dinal  intervals. 

(2)  Step  1  is  repeated  for  each  span  station  at  which  pressure  coeffi¬ 
cients  are  calculated. 

(3)  The  resulting  forces  per  station  are  then  plotted  against  their  par¬ 
ticular  span  station  to  again  calculate  the  area  under  the  curve  which  results 
in  the  total  induced  force.  This  procedure  is  followed  for  each  rectangular 
segment.  The  resultant  pressure  coefficient  and  induced  lift  of  each  segment 

is  then  summed  for  all  segments  to  obtain  the  induced  lift  for  the  total  planform. 
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The  accuracy  of  this  procedure  is  dependent  upon  the  integration  interval 
chosen.  This  is  especially  true  in  close  proximity  to  the  jet  exit  when  the 
surface  pressure  variations  are  quite  large  for  only  small  changes  in  location. 
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Y 

NOTES: 


L.  Rectangular  segments  are  smaller  near  the  location  of  the  jet  than  those 
farther  away.  Some  consideration  must  be  given  in  determining  the  inte¬ 
gration  intervals. 

2.  Rectangular  segment  lengths  are  sized  for  accuracy  in  modelling  the  plan- 
form  contour. 

3.  For  applicable  configurations,  only  half  of  configuration  panelled,  with 
results  multiplied  by  two  for  symmetry. 

4.  Jet  center  is  considered  to  be  the  origin  with  the  central  most  panel (s) 
offset  from  the  origin  to  avoid  a  zero  value  for  Y/D. 

5.  All  lengths  are  non-dimensionalized  by  the  jet  diameter.  Accordingly, 

sli  input  data  for  any  configuration  should  be  similarly  non— dimens ionalized . 


FIGURE  A-l.  Schematic  of  Configuration  Panelling  Procedure 
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Computer  Program 

An  interactive  computer  program  was  developed  to  facilitate  the  calcula¬ 
tion  of  the  jet  induced  pressure  coefficient  over  any  conf iguration  planform. 
The  program  contains  equations  (1)  to  (4)  with  which  to  calculate  the  pressure 
cofficients  for  any  given  point,  the  methodology  with  which  to  panel  the  con¬ 
figuration  planform,  and  the  integration  routine  to  calculate  the  induced  lift 
from  the  pressure  distriixition.  Programing  these  equations  has  also  enabled 
calculation  of  the  pitching  moment  resulting  from  the  planform  pressure  distri 
bution  according  to  equation  (A-l ) 

n 

AM  =  S  /atA  Xi 

TOi  1=1  \T  h  ^  (A-D 

where 


(AL\  ■  induced  lift  of  a  particular  panel 
T  }  i 

=  the  X  location  of  the  panel  center  of  pressure 


n  =  number  of  panels  used  to  simulate  the  planform 


The  planform  center  of  pressure  is  then  calculated  according  to  equation  (A-2 


(A- 2) 


As  output,  the  program  provides  the  pressure  coefficient,  induced  lift, 
pitching  moment,  and  center  of  pressure  for  each  panel  used  to  model  the 
planform  and  for  the  entire  planform.  A  listing  of  the  program  along  with 
a  sample  input  and  resulting  output  is  contained  at  the  end  of  this  appendix. 
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The  input  required  by  the  computer  program  consists  of  data  required  to 
panel  the  configuration  planform,  the  jet  location,  and  the  equivalent  velocity 
ratio.  This  input  is  tabulated  below. 


COMPUTER  INPUT 


ITTL  Case  Title  (Up  to  80  Characters) 

NCOMP  Number  of  rectangular  segments 

NVE  Number  of  velocity  ratios 

VE  Velocity  ratio 

DELXY  Integration  interval  into  which  rectangular  segments  are  divided 
(interval  for  both  x  and  y  directions;  different  values  can  be 
input  for  the  various  segments) 

X0D1  Minimum  x-coordinate  of  rectangular  segment 

X0D2  Maximum  x-coordinate  of  rectangular  segment 

Y0D1  Minimum  Y-coordinate  of  rectangular  segment 

Y0D2  Maximum  Y-coordinate  of  rectangular  segment 

X0D1,  X0D2,  Y0D1 ,  Y0D2  must  be  repeated  for  each  rectangular  segment. 
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